Qualitative analysis of dissipative cosmologies 



A. Di Prisco* 
Escuela de Fisica, Facultad de Ciencias, 
Universidad Central de Venezuela, Caracas, Venezuela 

L. Herrera^ 
Area de Fisica Teorica, Facultad de Ciencias, 
Universidad de Salamanca, 37008 Salamanca, Spain 

J. Ibanez-f 
Departamento de Fisica Teorica, 
Universidad del Pais Vasco, Apdo. 644, 48080 Bilbao, Spam 
(February 7, 2008) 

The evolution of an homogeneous and isotropic dissipative fluid is analyzed using dynamical 
systems techniques. The dissipation is driven by bulk viscous pressure and the truncated Israel- 
Stewart theory is used. Although almost all solutions inflate, we show that only few of them can be 
considered as physical solutions since the dominant energy condition is not satisfied. 
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I. INTRODUCTION 



There has been a renewed interest in the study of dissipative cosmologies related to the existence of inflationary 
solutions. As a matter of fact, dissipative processes are thought to be present in any realistic theory of the evolution 
of the universe. The simplest model of dissipative fluid is that which assumes the existence of bulk viscosity only. At 
a phenomenological level, bulk viscosity may be associated with particle production M, ||. Also, qualitatively, bulk 
viscosity may be understood as the macroscopic expression of microscopic frictional effects that appear in mixtures 
H. The easiest way to include bulk viscosity effects is through the Eckart's theory j|, which assumes that the bulk 
viscous pressure is proportional to the expansion. Several authors have used that theory to investigate the effects of 
viscosity on cosmological models || . 

As it is well known, however, Eckart's theory suffers important drawbacks || and a more complete theory must 
be used. One of the proposed theories is the Israel-Stewart theory Jjj which complies with causality and stability. 
The non-linear terms of that theory are often neglected giving place to the so called truncated theory. Although the 
truncated theory may lead to a different behaviour as compared with the full theory, we shall for simplicity base our 
analysis on the truncated version of the Israel-Stewart approach. This truncated teory was firstly used by Belinskii 
ct al. Ji| in a cosmological context and after that by many authors ( see Q and references therein). The flat FRW 
models were studied in ref. [jl0| and has been extended recently to spatially curved solutions [jllj . In these last works 
the system of field equations was reduced to a second order differential equation and the stability of the stationary 
solutions was investigated by using a Lyapunov function. A different and more powerful approach was given by Coley 
and van den Hoogen ]T^ | who used dynamical system techniques to analyze the asymptotic behaviour of the governing 
field equations. They used dimensionless equations of state that make the equation for the expansion to decouple 
from the system, obtaining a two dimensional system of differential equations. 

Since these models were originally studied, it was soon realized that the type of equations of state used in them 
was determinant in verifying whether the models underwent inflation. It is therefore important to study the effect 
that other reasonable equations of state will produce in the qualitative behavior of the solutions. 

In this paper we analyze the asymptotic behaviour of isotropic and spatially homogeneous models with a fluid with 
bulk viscosity It has been pointed out J(| that although bulk viscosity is almost vanishing, in both the ultrarelativistic 
and the Newtonian regimes, its dynamical effects can not be neglected. On the other hand the evolution equations 
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given by this type of matter can be used to model another kind of sources (e.g. the string dominated universe described 
by Turok @ or even a scalar field) The equations of state we assume are those introduced by Belinskii et al. ||. 
Unlike the case studied in ref . , now the expansion does not decouple and we should deal with the complete three 
dimensional system, except for a particular value of the parameters. It is important to stress that the choice of the 
equations of state could dramatically change the behaviour of the model and, therefore, it is worthwhile to study the 
implications of different reasonable equations of state. 



II. FIELD EQUATIONS 



Let us consider an homogeneous and isotropic spacetime with an imperfect fluid moving orthogonally to the hyper- 
surfaces of constant curvature. The source is a fluid with bulk viscosity: 



T ab = pu a u b + (p + a) (g a b + u a u ) , 



(1) 



where a is the viscous pressure p is the thermodynamic pressure and p is the energy density. The Einstein field 
equations and the conservation equations are written as: 



p = -3H(p + p + a) 
H = -H 2 - -(p + 3p + 3a) 
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3H 2 = p-^, (A = ±1,0), 



(2) 
(3) 

(4) 



where a is the scale factor and H = — is the Hubble parameter. We assume an evolution equation for the viscous 
pressure a given by the so called truncated Israel- Stewart theory (see |l4j and references therein): 

a + t& = -3£H. (5) 

where £ is the coefficient of bulk viscosity and r is the relaxation time (£ > 0, t > 0). Although eq. (5) (strictly speaking) 
is valid only on the assumption that the fluid is close to equilibrium (|er| <C p) we will assume, for simplicity, that it 
describes the evolution of the viscous pressure even for regimes far from equilibrium. A more complete description 
would require taking into account the full transport equation of the Israel-Stewart theory. We assume a barotropic 
equation of state p = (7 — l)p, < 7 < 2 and the following relations: 



£ = ap m , t = —, to, a > 
P 

which were used by Belinskii et al. S. 

In order to study the system of equations we define adimcnsional variables and a new time: 



(6) 



P 

3# 2 ' 



17 

H 2 ' 



H(t)dt = dr. 



Now the equations (2), (3) and (5) reduce to 

Cl' = (0- l)[fi(37-2) +£] 



H' 



-90 



-H 



1 



37 



Q 

- 2 



-fi 



1. 



2- (37-2)fi-E 



(7) 

(8) 
(9) 

(10) 



2 2 

where ' means derivative with respect to r. In addition to the above equations we have the Friedmann equation (4) 



(ii) 



We will assume H positive in the open and flat models, so the system (8)-(10) is well defined. This is not true in 
the closed case, for which H could become zero. So in this case the system should be understood as describing the 
evolution of closed models at early times when the expansion is positive. 
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Before studying the qualitative behavior of the system it is important to consider the energy conditions. The strong 
energy condition (sec) implies that p + 3p e ff > 0, where p e ff = P + 0". In terms of the new variables this condition 
is written as 

(3 7 - 2)fi + £ > 0. (12) 

The violation of this condition implies that the scale factor accelerates and the solution inflates. The dominant energy 
condition (dec) requires p + p e // > or: 

37« + S>0. (13) 

Matter not verifying this condition is considered unphysical, futhermore, it seems that the violation of this condition 
could lead to the violation of the generalized second law of thermodynamics O . 



III. QUALITATIVE BEHAVIOUR 

The line fl = 1 is an invariant space and divides the phase space of the system into three invariant sets: fl > 1 
corresponding to closed models, O < 1 corresponding to open models and Ct = 1 to flat FRW solutions. Besides these 
invariant spaces, H = is another invariant set that splits again the phase space into another three invariant sets. 
However we do not consider the H < space for we assume that in early times H > 0. 

When to = 1/2 the equation (10) decouples from the two other equations. Therefore the system becomes a two 
dimensional system which can be studied much more easily. To investigate the qualitative features of the system we 
will start with this case and then with the complete system. 
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In this case the system reduces to a two dimensional system of the form: 

fl' = (Q- l)[fi(37-2) + E] 



£' = -90 - £ 



'30 



- 2 - (37 - 2)0 - £ 



(14) 
(15) 



There are four equilibrium points: 



= 1, £ = £* = - 




P 3 : O = 0, £ = 0, 



(16) 



/ a 67 + 5 
V7I37-2 



,£ = -(37-2)0, 7 >-. 



The stability of these points is given by the sign of the eigenvalues of the matrix of the linearized system around each 
point. The point P3 has two eigenvalues of different signs, so it is a saddle point. The two eigenvalues of the point 
P + are positive and, therefore, it is a source point. The point P~ has its two eigenvalues negative provided a > cxq, 
where 



ao 



_ V3(3~f - 2) 
67 + 5 : 



(17) 



being, therefore, an attractor. When a < a this point becomes a saddle. The behaviour of points P4 and P - are 
opposite: when a > ao, Pa is a saddle but is an attractor when a < ao. Let us note that since a > the condition 
a < ao implies 7 > 2/3 and £lo < 1. Therefore in this case the only attractor is the point P4. When a — ao both 
points, P4 and P~ coincide. 

The exact solutions corresponding to each points are the following: P3 is the Milne universe for which a(t) = t. At 
the points P ± we have 



3 



a(t)=a t H \ p(t) = ^-, a(t) = ^-^-, H = - * g± . (18) 

At P + we have Ho < 1 so the solution is a FRW non inflating model. At P~ the solution inflates when a > ao 
(Ho > 1) which is equivalent to say that the solution does not verify the strong energy condition (12). The solution 
at P 4 is: 

/67 + 5a\ 2 l . . 37 — 2 . . 

O(t) = 0o*, p(t) = 3 ^g^— ^ 7y= J ^, <K*) = ^— P, (19) 

and the deceleration parameter g = 0. 

It is important to check whether the dominat energy condition (13) is satisfied by the solutions. This condition is 
always satisfied at the points P + and P4. It is easy to see that when the point P~ is an attractor (a > ao) then eq. 
(13) is verified as long as 7 > a/y/3. When P~ is a saddle, a < ao means that 7 > a/y/3 and the DEC is verified at 
this point. In Table I the main results concerning points P4 and P~ are summarized. 

To conclude with the description of these solutions we will compactify the phase space in order to analyse the 
behaviour at the infinity. To do this we introduce polar coordinates: SI = rcos0 and £ = r sin0. To make r finite we 
define 



1 

and a new time r: 

dr 



(20) 



a-r l ~ r - (21) 



With these definitions the system (18)-(19) can be written as 

dr 
~dT 



( 1 - r) ( r 2 [(37 - 2) cos 9 + sin 01 - - y/3r(l - r) cos 9 sin 2 9 
I a 

+(1 - r) [2r - 3~/r cos 2 - lOr cos 9 sin 9} } 



rlf) 

— = (1 - r) [-9 cos 2 + sin 2 9 + 37 cos sin 01 



- - J3r(l - r) cos sin cos 0. (22) 
a 

We deduce from the above equations that the circle r = 1 consists of equilibrium points. Close to the circle we get 

^ ~ (1 - r)r 2 [(37 - 2) cos + sin 01. (23) 
or 

When 37O + S > the above derivative is positive and, then, the points in the circle are attractors, while when 
3757 + S < the derivative is negative and the points are sources. 

In Fig.l we have plotted integral curves of the system (22) for different values of the parameters. The dashed curve 
represents the invariant space SI = 1 . Points to the left of each vertical lines violate respectivaly each of the energy 
conditions. Curves that start or end at 9 — ±7r/2 are unphysical since Q becomes negative and that means that at 
some time of their evolution the energy density is negative. 

In Fig. 1(a) the orbits starting at P + and going to the point P~, which is the only attractor for these values of the 
parameters, enter into the region where the strong energy condition is violated, so the solutions start inflating. Since 
P~ violates the DEC, the solutions enter into a region where this condition is violated being, therefore, unphysical. 
So, we see that in this case the only physical solutions are those that start at P + and go to infinity. These physical 
solutions do not inflate. In Fig. 1(b) we have plotted orbits corresponding to values of the parameters such that the 
point P4 is a saddle point. The main features in this case are similar to that showed in Fig. 1(a). The main difference 
is that now, since the attractor P~ verifies the DEC, then there are physical solutions ending at this point. Finally in 
Fig. 1(c) we have plotted orbits corresponding to the case for which the attractor is the point P4. What is remarkable 
in this case is that there are solutions that after undergoing inflation their final state is a universe with Q < 1 showing 
that the presence of a state of inflation in the universe does not imply necessarily that J! = 1 jlj| . 

Let us notice from these figures that if we assume a solution starting close to the equilibrium regime, i.e. | cr] <C p 
or |S| ~ then, at the end of the evolution, |E| is of the same order than SI which is out of the range of validity of 
the theory. This is the so called "runaway" solution Q . There are however solutions that fulfil this requirement. If 
we take a<Cl then E~ ~ — 3\/3a + 0(a 2 ) and, in this case, we can choose the parameters in such a way that P~ is 
an attractor that verifies DEC and the physical pressure p + a of the solution is positive. So, solutions starting with 
I XJ j ~ evolve towards P~, being close to equilibrium 
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B. m/ 1/2 



In this case the evolution is described by the three-dimensional system (8)-(10). Although it is more difficult to 
visualize the evolution, the behaviour is similar to that described in the previous section. To further simplify the 
equations we define a new variable 



h = H 1 ~ 2m . 



The system (8)-(10) is cast into a simpler form 



Q' = (fi-l)[fi(37-2) + E] 
(SIT) 1 -" 1 



£' = -9Q - s 
ft' = -(1 - 2m)ft 1 + 
The equilibrium points of the above system are: 



ft - 2 - (37 - 2)fi - £ 



3 7 



2 2 



(24) 

(25) 
(26) 

(27) 



P ± :0 = 1, ft = 0, S ± = -(- 7 ± Vl 2 + 4) 
P 3 : n = 0, ft = 0, E = 0, 

P 4 :0 = lft= , £ = -3 7 . 

7 

In addition to these points there is another one, P5, which is valid only when m > 1/2: 

P 5 : = 1, P = (ft = 00), E = 0. 



(28) 



(29) 



Let us note that when to > 1/2 the equilibrium points, except P4 and P5, have P = 00. 

As before, by looking at the eigenvalues of the linearized system we obtain the stability of these points. The points 
P ± have eigenvalues 



3 7 - 2 + E= 



3 7 + 2S ± , 



1 - 2to 



(37 + S±). 



(30) 



For P + the two first eigenvalues are positive and the third is negative when m < 1/2 and positive when to > 1/2. 
Thus, this point is a saddle point if to < 1/2 and a source if to > 1/2. For P~ the two first eigenvalues are negative 
and the third is positive when m < 1/2 and negative when m > 1/2. Thus, this point is a saddle point if m < 1/2 
but is an attractor if to > 1/2. 

The eigenvalues of the point P3 are: 



-(4-37 + 3^7 + 4), -(4-37-3VA7T4), 



-(1 - 2m). 



(31) 



The first eigenvalue is positive and the second negative. The third is either negative or positive depending on the 
value of to. Thus this point is a saddle point. Finally P4 eigenvalues are 



-2, 



2 7 



[-( 7 2 + l)± V / (7 2 + 1) 2 -2(1-2to)]. 



(32) 



If m < 1/2 the three above eigenvalues are negative and, then, P4 is an attractor. But if to > 1/2 one of the 
eigenvalues is positive and P4 turns to be a saddle point. 

To analize the stability of the point P5 we change to spherical coordinates: O = rsin#cos0, £ = r sin 9 sin <j>, h = 
rcosO and we compactify the phase space by using the coordinate r given by (20) and the time r defined by (21). 
Close to the point P5, (r <~ 1, sin 9 ~ 1 — r, (f> = 0) we obtain: 



dr 1 

-~--(l-2TO)3 7 (l-r) 2 . 



(33) 



The above expression is positive when to > 1/2 and therefore this point is an attractor. 
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The solution described by the point P4 is given by 



/om \ 1/(1-2to) 

a{t) = a e Hot , p = 3P 2 , a = -3 7 P 2 , H = f — J . (34) 

This de Sitter solution was already obtained by Barrow [Q. Since the rest of equilibrium points have H — (or 
infinity depending on the value of m) it is not possible to get the corresponding exact solutions (except for the point 
P3 that corresponds to the Milne solution) We can, however, obtain approximate solutions for some of the points. 
When to > 1/2 the point P + is a source. By linearizing the system (25)-(27) around this point we get 

t^-oo, H~H a e- XT , A=i(3 7 + £+), (35) 

where Hq is a constant of integration. Integrating the other two variables and using (11) to change from r to t we 
obtain: 

t-0, p~^p a-^, «(f)~*V\ (36 ) 
To see the behavior close to the point P5 we can integrate eq.(27) and after a straightforward calculation we obtain 

*->«>, H ~^- ( 37 ) 

This behaviour in the vicinity of the point P5 was already obtained by Belinskii et al. S and by Barrow 

As to the energy conditions, it is easy to see that both are satisfied by the points P^ and P3 and are violated by 

the point P~ . The DEC is satisfied by P4 and P5 but SEC is not satisfied by P4 and is satisfied by P5 when 7 > 2/3. 
The behavior described above is summarized in Table II. 

The compactification of the phase space in the case to 7^ 1/2 is much more difficult than in the m = 1/2 case, 
making almost impossible to get a complete description of the system. However we can obtain the main features 
concerning the qualitative evolution of solutions from the information on the equilibrium points and the numerical 
integration of the system. 

In Fig. 2 we have plotted the projections of the orbits, obtained by numerical integration of the system (25)-(27), 
onto the plane (f2, £) for different values of the parameters. The straight lines indicate whether the solutions fulfil 
each of the energy conditions. Points below each line do not satisfy the corresponding energy condition. As in the 
former case there are solutions that enter to or come from the f2 < region. All these solutions are unphysical. 

Fig. 2 (a) corresponds to the case where the only attractor is the point P4. Besides the solutions that end at 
infinity, there are solutions approaching point P + , which is a saddle, then entering into an inflationary regime and as 
approaching the point P _ , which is a saddle as well, they become unphysical since the dominat energy condition is 
violated. The values of the parameter in Fig. 2(b) are such that there are two attractors: points P~ and P5. Since in 
this case P~ does not satisfy the DEC, the solutions going to this point become unphysical. However, solutions ending 
at P5 could be physical. These last solutions could be of physical interest for they are the only solutions that could 
fulfil the requirement that the viscous pressure be much less than the equilibrium pressure during their evolution. 
Starting at points for which £ ~ they evolve towards solutions with £ = 0. 

IV. CONCLUSIONS 

In this paper we have analyzed the qualitative behaviour of a particular type of dissipative fluids in FRW universes 



using the truncated Israel-Stewart theory. There have been a previous |10| analysis of the same system and with the 
same equations of state but only for the zero curvature FRW universes that was thereafter extended to the general case 
jLl| . In those analysis the system of equations was reduced to a second order differential equation whose stationary 
solutions were studied by using a Liapunov function. That method does not discriminate, however, between different 
equilibrium points. We have seen in this paper that when to 7^ 1/2 almost all the equilibrium points have H = and 
all of them were reduced to only one stationary solution in the mentioned papers. 

An important feature that emerges from our analysis is that for most of the equilibrium points either the dominat 
energy condition or the positivity of the energy density are not satisfied. Solutions evolving towards these equilibrium 
points, therefore, should be ruled out as unphysical. That means that in order to obtain physical solutions we have 
to be very careful with both the initial conditions and the values of the parameters. With respect to inflation, we 
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see that, generically, the presence of bulk viscosity makes the solutions enter into an inflationary regime. Only few 
solutions do not inflate. 

If we consider the fluid interpretation of our system of equations, we have seen that almost all the solutions evolve 
in such a way that they go outside the regime in which the theory is valid (i.e. when |cr| Cp). The only exceptions are 
those solutions which evolve towards the point P5 when m > 1/2 which describes an inflating universe when 7 < 2/3 
or solutions evolving towards P~ with m = 1/2 and a very small. 

As mentioned before we have used the truncated equation (5), instead of the more general transport equation 
derived from the Israel-Stewart theory. As it has been shown p4| , this latter equation may be written as 



-3^H 



1 

7cg \P 



(38) 



with 



and 



c 2 - * 

c b — 



{p + P)l 



1 + 37Tff 



; £* 



1 + 3jtH 



(39) 



(40) 



Therefore, close to equilibrium [a << p), equation (38) leads to the truncated equation (5) with reduced relaxation 
time (t*) and bulk viscosity (£*) 



T*<7 + cr = —3£*H. 



(41) 



The ammount of reduction depending on the size of £ relative to H. Obviously if tH << 1, then t» w t and « £. 
However if tH is close to 1, the reduction could be significant. On the other hand, qualitative changes in the results, 
outside the quasi-equilibrium regime, are difficult to forecast without considering specific models. 

It is important to realize that the system of equations studied in this paper, not only describes the evolution of a 
dissipative fluid with bulk viscosity, but also the evolution of a different kind of matter. It is well konwn, for instance, 
that when m = 2/3 the above equations give the evolution of a string dominated universe jl3j . We shall see that this 
equivalence can be extended to a scalar field as well. 

Indeed, a scalar field in both FRW and Bianchi type solutions is equivalent to a perfect fluid with p and p defined 
as: 



p = \& + v(4>), 



(42) 



where V is the scalar field potential. We can substitute in the above expression p by p + a. If at the same time 
we consider a barotropic equation of state between the energy density and the equilibrium pressure of the form 
p = (7 — l)p then we obtain that the viscous pressure is given by: 



' 7 



-7V. 



(43) 



Taking derivatives with respect time and substituting into the Klein-Gordon equation of the scalar field 

4> + 3H<j) + = 
d(p 

we obtain that a satisfies an equation like (5) with r and £ given by 



— f = — ( V — — 

6iT ^ ~ 3H V + 3H dd> ' 



(44) 



(45) 



Now, in dissipative systems there exists a peculiar state called "critical point" , which corresponds to a specific value 
of a combination of thermodynamic variables. It has been shown |l6|, |l^] (and references therein) that a dissipative 
system in the critical point immediately (i.e. on a time scale of the order of the relaxation time) after its departure 
from equilibrium, behaves as if the effective inertial mass of any fluid element vanishes. 
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2(p + p). (46) 



In the case of pure bulk viscosity the critical point is characterized by 

£ 

T 

In this case however the critical point is forbidden by causality and stability requirements, demanding 



<(p+p)(l- c 3) (47) 

where c s is the speed of sound (see fl6f|). 

Nevertheless in order to see what the critical point may imply in terms of a scalar field, let us assume (^6|). Then, 
combining with (|4^), we obtain after an elementary algebra 

1 dV UK) 

3H^~T (48) 

and feeding back into ([13) 

+ 3H0(l + !) =0. (49) 

A simple integration of ( fl9"| ) gives 

^ e -/ 3ff(1+7/2)dt . (50) 

Let us consider two examples: 

• de Sitter case, H =constant. 
In this well known case, one gets from (Eih 



-3H(l+7/2)i 



and therefore 



Then (ph yields 



(51) 



(52) 



V~<j) 2 . (53) 



This kind of potential is one of the most favoured by particle theorists, since it describes renormalizable particle 
theory @. 

• H~ ^- 

37* 

The asymptotic behaviour close to the point P5 described above (see eq.(37)), yields, using ( |5C| ) 

~ ^ /2) (54) 

and therefore the corresponding potential is 

y-0 7+2 . (55) 

Thus, the critical point condition (fl6|), leads, for an scalar field interpretation of the source term, to polinomial 
potentials (for the two examples considered). 

It is evident, therefore, that the constraints imposed by the fluid interpretation can be droped if we consider the 
energy momentum tensor to correspond to another kind of matter (different from the fluid with bulk viscosity). In 
such a case the only warnings to bear in mind are those coming from the dominat energy condition requirement. 
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FIG. 1. Plots of the integral curves of the sytem when m = 1/2. Vertical axis represents the coordinate r an the horizontal 
axis represents the 9 coordinate. They range from to 1 and from —n/2 to n/2 respectively. Dashed line represents the curve 
Q = 1. Points to the right of the SEC line verify the strong energy condition and points left don't. The same for DEC line, 
corresponding to the dominat energy condition, (a) In this case 7 = 1/3 and a = 1, (a > ao, 7 < a/\/3). (b) Now 7 = 2 and 
a = 1, (a > a ) and (c) 7 = 2 and a = 0.2, (a < a ) 



FIG. 2. Projections of the integral curves of the system when m 7^ 1/2. Points below the solid diagonal line does not verify 
the strong energy condition, while those below the dashed diagonal line does not verify the dominat energy condition, (a) In 
this case 7 = 1/3, m = 1/4 and a = 1. (b) Now 7 = 1/3, m = 0.7 and a = 1. 



TABLE I. Qualitative properties of points P and Pi 





Pi 




P~ 




DEC 












7 > q/V3 


Yes 


a > ao 


7 < 2/3 




Attractor 




7 < a/y/H 


No 












7 > a/V3 


Yes 


a > eta 


7>2/3 


Saddle 


Attractor 




7 < a/s/3 


No 


a < ao 


7>2/3 


Attractor 


Saddle 






Yes 
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TABLE II. Qualitative properties of points equilibrium points 





P+ 


P~ 


P 3 


Pa 


P 5 


m < 1/2 


Saddle 


Saddle 


Saddle 


Attractor 




m > 1/2 


Source 


Attractor 


Saddle 


Saddle 


Attractor 


DEC 


Yes 


No 




Yes 


Yes 


SEC 


Yes 


No 




No 


if 7 > 2/3, Yes 
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